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ABSTRACT 


The dynamics and stability of large orbiting flexible beams, and plat- 
forms and dish type structures oriented along the local horizontal are 
treated both analytically and numerically. It is assumed that such 
structures could be gravitationally stabilized by attaching a rigid light- 
weight dumbbell at the center of mass by a spring loaded hinge which also 
could provide viscous damping. For the beam it is seen that the small ampli- 
tude inplane pitch motion, dumbbell librational motion, and the anti-symmetric 
elastic modes are all coupled. Tlie three dimensional equations of motion 
for a circular flat plate and shallow spherical shell in orbit with a two- 
degree-of freedom glmballed dumbbell are also developed and show that only 
those elastic modes described by a single nodal diameter line are influenced 
by the dumbbell motion. Further, in the case of shallow spherical shells 
the pitch and the axi-symmetric modes are seen to be weakly coupled in the 
linear range. With the shell’s symmetry axis following the local vertical, 
the structure undergoes a static deformation under the influence of gra- 
vity and inertia. Stability criteria are developed for all the examples 
and a sensitivity study of the system response characteristics to the 
key system parameters is carried out. 
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CHAPTER 1 


INTRODUCTION 

This represents the second part of the 1979-80 final report and con- 
centrates on the modelling, dynamicil', and stability of large, flexible 
earth pointing orbiting stmctuteSr* First, a brief review of the general 
formulation (continuum model) for any large flexible space system in 
orbit, as previously developed*^* is summarized in Chapter 2. 

The paper to be presented at a forthcoming conference forms the basis 
of Chapter 3: '*0n the Dynamics and Control of Large Orbiting Flexible 

Beams and Platforms Oriented along the Local Horizontal,” XXXIst Congress 
of the International Astronautical Federation, Sept. 21-28, 1980, lAF- 
80-E-230. This section treats the uncontrolled dynamics of a large thin 
flexible beam in orbit with emphasis placed on the motion about tlie local 
horizontal orientation instead of the J^ocal vertical. The use of a 
gimballed connected dumbbell is proposed to provide the correct composite 
moment of inertia ratio required for gravitational stabilization and to 
also offer a restoring torque (and possibly damping torque) due to the 
spring-glmball assembly. A further extension of this concept considers 
the use of a two-degreei-of-freedom gimballed dumbbell to aid in the sta^ 
bilization of a large flexible plate (platform) in orbit about the nominal 
orientation In the Jovial horizontal (tangent) plane 

In Chapter 4 the treatment of Chapter 3 is extended to consider the 
dynamics and stability of a large flexible shallow shell structure in 
orbit. As in the previous formulation for the beam and plate, a sta- 
bilizing dumbbell with two degrees of freedom is also incorporated to 
offer gravitational stabilization characteristics. 

In Chapter 5, general concluding comments and recommendations for future 
work to be initiated in the 1980-81 grant year, in accordance with a 
recent proposal^ to NASA are discussed. 


For references cited in this report, see list of references 
provided after each chapter. 
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aiAPTER 2 


EQUATIONS OF MOTION OF AN ARBITRARY FLEXIBLE BODY 
IN ORBIT 

A brief review of the basic equations describing the attitude motion 
of an arbitrary flexible body in orbit is presented here. Details 
regarding the derivation of these equations can be found in ref. 1. 

In general, any arbitrary motion of a flexible body in space, ?ibout its 
center of mass, consists of a rigid body rotation about an Instantaneous 
axis of rotation through the center of mass and elastic deformations. 

The rigid body rotation can be expressed by a sequence of three successive 
Euler angle rotations about the body fixed ax^^a. Furthermore, if we 
assume small elastic deformations of the body, the elastic displace- 
ments can be expressed as the superposition of an infinite number of 
structural mode shapes weighted by the time dependent amplitude functions. 
Expressed symbolically, this appears as, 

q (rQ,t) - A^(t) <I> (rQ> (2.1) 

where, Icq » position vector of a generic point in the body in the 
undeformed state ^ 

q “ elastic displacement vector at r^ 

$ » vector of mode shape functions. 

Thus, the basic equations which describe the general motion of an arbitrary 
flexible body in space about Its center of mass consist of (n+3) equations: 
three equations describing rigid body rotational motion; and n-modal equ- 
tlons which describe the elastic motion of the body. These equations 
as developed in ref. 1 appear as follows; 

Equations of rigid body rotation; 

R + S Q E D^”^ » + E + C (2.2) 

where, R = inertia torques due to rigid body motion 

m inertia torques due to elastic motion 

j,p(n) » torques due to center of mass shift effects 
n (zero for uncons trained structures) 


Ref. 1; Bainum, P.M. , Kumar, V.K. , James, P.K. , "The Dynamics and 
Control of Large Flexible Space Structures," NASA CR-156976 Report, 
May 1978. 
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■ gravitational torques due to rigid body motion 

" gravitational torques due to elastic motion 
C " other external disturbance and control torques 


Equations of elastic motion: 
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The term R in eq. (2.2) contains the terms from the classical Euler 
equations for rigid body rotations which are obtained by assuming the 
^ody to be rigid. In the body fixed principal axes reference frame, 
R is given by, 
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where, (J^,J ,J^) « principal moments of inertia of the body in 
^ the undeforraed state 

(a)^,U3 • components of the body angular velocity vector, 

^ w, in the body principal axes frame 

A i\ A 

l,j ,k, " unit vectors along the body principal axes 


The vector expressions for the other remaining terms ixi eqns. (2.2) and 
(2.3) are as follows: 
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where, a ■ Inertial acceleration of the center of mass of 
the body 

fQ ■ gravitational force per unit mass at the center 
of mass of the body 

r » instantaneous position vector of a generic mass 
point in the body 

e * external disturbance and control forces per unit 
mass of the body 

M ■ a matrix operator dependent on the Euler angles 

Tlie mode shape vectors, (n » 1,2,...) are orthogonal to each 

other, i.e, 

/ (2.15) 

V n mn ' 

In addition, for unconstrained structures in space, rigid body translation 
and rotational modes should be orthogonal to (n * 1,2,...) 

/ - 0 (2.16) 

V 

/ rQX<I>^"^dm » 0 (2.17) 

V 

Several specific applications of flexible spacecraft in orbit considered in 
the following chapters utilize the equations presented in this chapter as 
the basis of the foraulatlon of their models. 
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Chapter 3 


ON THE DYNAMICS OF LARGE ORBITING FLEXIBLE BEAMS AND 
PUTFORMS ORIENTED ALONG THE LOCAL HORIZONTAL 


ABSTRACT ‘ 

The dynamics and stability of large orbiting flexible beams and platforms oriented 
along the local horizontal are treated both^analytically and numerically. It is 
assumed that such structures could be gravitationally stabilized by attaching a 
rigid lightweight dumbbe.ll at the center of mass by a spring loaded hinge which 
also could provide viscous damping. For the bean it is seen that the small ampli- 
tude inplane pitch motion, dumbbell lib rational notion, and the anti-sjrmmetric 
elastic modes are all coupled. Tlie three dimensional equations of motion for a 
circular flat plate in orbit with a two-degree-of-freedora giraballed dumbbell are 
also developed and show that only those elastic nodes described by a single nodal . 
diameter line. are influenced by the dumbbell motion. Stability criteria are de- 
veloped for both examples and a parametric study of the least damped mode charac- 
teristics together with numerically simulated transient responses are carried out. 


KEYWORDS 

I 

Flexible spacecraft dynamics; stability; stabilizing gimballed dumbbell booms; 
large space structures, 

t 

k 

INTRODUCTION 

Proposed future applications of large space structures include: space based power 

generation and transmission (to earth); communications; earth resource observation 
missions; and electronic mail systems. 
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Through the use of such systems, one can see an example of the applications of 
space developments to solve some of the ptoblems of mankind such as in education, 
economy, and energy. The applications described here all require that the largest 
surface (length) of the system be nominally oriented along the orbital tangent or 
normal to the local vertical. To gain insight into the types of problems in- 
volved with inherently complex space structures, two simple examples of large flex- 
ible space systems are addressed in this paper. It should be recognized that beams 
and plate elements would be two of the most fundamental structural elements in any 
large apace structural system. 

Previously the equations of ^motion of a general arbitrary flexible spacecraft in 
orbit have been developed. The elastic displacement at any arbitrary point in 
the structure was assumed to result from a superposition of the different flexural 
modes. The various terms in the general vector equations of motion were expanded 
in terms of the spatially dependent modal shape functions and frequencies As 
a specific example, the dyna^cs and stability of a long, flexible beam constrained 
to move only in the orbital plane was considered, with the principal emphasis 
placed on the motion about the nominal earth pointing (local vertical) orienta- 
tion. ^ It was observed that for small amplitude pitch and flexural oscillations, 
the pitch motion was not affected by the elastic modes, and that the elastic mo- 
tion was coupled to the pitch motion and described by sets of Mathieu equations. 
The possibility of parametric instability at very low natural elastic frequencies 
was demonstrated. 

The present paper extends the work of Ref. 2 to consider: (1) the motion and 

stability of the beam about a nominal local horizontal orientation; and (2) the 
dynamics of free-free homogeneous platforms in orbit with emphasis placed on a 
circular plate structure. 


DEVELOPMENT OF EQUATIONS OF MOTION-BEAMS ORIENTED ALONG ’ 

THE LOCAL HORIZONTAL 

A Uniform Beam in Orbit with its Axis Nominally along the Local Horizontal 


Fig. 1 shows a long, thin flexible beam in orbit with its centroidal axis nomi- 
nally along the local horizontal. The following assumptions are made in de- 
riving the equations of motion: (a) the. beam is long and slender with uniform 

cross section and uniform distribution of mass and stiffness properties; (b) the 
center of mass of the beam follows a circular orbit; (c) all motions and deforma- 
tions of the beam are res-tricted to occur 'within the orbit plane; (d) there are 
no constraints on the beam's elastic motion; (e) longitudinal vibrations of the 
beam are negligible in comparison to the transverse vibrations. • 


Based on assumption (c) , the Euler angles representing the out-of-plane beam roll 
and yaw motions vanish and also the out-of-plane conponent, i|)^’^,of the mode shape 
vector is set to zero. In addition as a result of assumption (e) the longitudnal 
component, of the modal shape vector also vanishes. Furthermore, for uncon- 

strained structures and by virtue of the orthogonality conditions together with 
the other assumptions, it can be shown that both the inertia and gravitational 
coupling terms between the rigid body modes and flexible modes (E E 

gj^) appearing in the rotational equations of motion and the generic modal equa’-> 
tions (Eqs. (15) and (17) of Ref. 2) vanish. In addition, the terms associated 
with the shift in the center of mass (S and also vanish for the assumed 

unconstrained motion. As a result, the rotational and generic equations of motion, 
for this application, simplify to; 
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( 1 ) 
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* 3£]j2 j sinO ; A 
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Jy » beam pitch axis moment of inertia 


flexural modal amplitude 

M “ n*-'^ modal mass 
n 




Cy = external torque about the pitch axis 


E = effect of external forces on the n 

n 


th 


mode 


for small amplitude pitch oscillations of the beam with respect to the local hori- 
zontal, sin0T'0, cos0=l, and Eqs.- (1) and (2) simplify to. 


0"-39 - Cy/Jyaj2 

®n (9'-l)2-2]e^ » -E^/M^w^i 


(3) 

(4) 


where 


e ~ A /a (I ~ undeformed beam length) 
n n 

(■ d/dt, where t = m^t-. 
t “ time 


It is well known that in the absence of external control torques, the nominal 
local horizontal orientation of the bcjam represents an unstable motion which is 
reflected by Eq. (3). In order to overcome the destabilizing effect of the gra- 
vity-gradient torque on the beam, one can either apply active control torques, 
or adjust the moment of inertia distribution of the system such that the gravity- 
gradient torque now becomes stabilizing. The beam can be gravitationally sta- 
bilized by using a rigid lightweight dumbbell with proper moment of inertia 
(Fig. 1). In this work the dumbbell is assumed to be attached at the center of 
the beam by a spring loaded hinge with viscous rotational damping also assumed 
to be presents The proper moment of inertia of the dumbbell can be attained by 
selection of the tip masses which are assumed to be much larger than the mass of 
the dumbbell rod. 
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Gcavltatlonally Stabilized Flexible Beam in Orblt-^Axis Nominally along 
the Local Horizontal 


In addition to the assumptions made in the development of Eqs, (1) (4), it 

will also be assumed that the stabilizing dumbbell is rigid and that the mass of 
the connecting link is negligible compared with that of the tip masses. In 
addition to the restoring torque provided by the spring, it will also be assumed 
that a linear viscous damping force is also provided at the hinge. 


Tlius, the resulting restoring and dissipative torque at the hinge can be ropre 
seated as, 


k(a-0-9) + c(a-o-O) 


(5) 


where 

k - torsional restoring spring constant at the hinge 
c *• viscous damping coefficient 

a « angle between the dumbbell axis and the local vertical (see Fig. 1) 

a " rotation angle of the beam normal'lit the hinge due to the elastic 
, deformation (Fig. 1) . 


For small elastic deformations (q ’ 'S A ), 
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(7) 


The modal component of the torque, Cy, can be obtained by replacing C by an 
equivalent force system consisting of two forces of equal magnitude but opposite 
in sign and separated by a small distance. In the limit as th'?> forces are moved 
closer, 
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n 



3z 


( 8 ) 


z«0 


The pitch and modal equations for the system shown in Fig. 1 are easily obtained 
by substituting Eqs. (5)-(8), into Eqs. (1) and (2), with the result, in dimen- 
sionless form, 

9”-3sin0 » k(a-0)+c(a'-e')-S(ke (9) 

n n n z 

e”+C (w^/w^)^-(0'-l)^-23e^ “ {k(a-0)+c(a' ~0 ’) 
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} /M 

z y ti 


( 10 ) 



vhere 


k - k/JyU)2 . c - c/JyW^ 

(n) 


,(n) 


d(^ 


X 




z«0 


The equatrion of tnoclon of the dumbbell is given by. 


a”+3sina " -kc, (ct-0) -cc- (a ' -0 ’ ) +JI (ke -1-ce ’ ) c. C 

JL ^ li «!• 


(n) 


( 11 ) 


where , 


'1 ■ 


“ pitch moment of inertia of the dumbbell. 

The following observations can be made from a study of Eqs. (9) -(11): (a) the 

pitch motion of the beam, the dumbbell motion (a) , and the elastic motion of the 
beam (e ) are all coupled to each other; (b) within the linear range the elastic 
modes for which » 0 (the symmetric modes) , arc completely independent of the 

pitch and dumbbell motions. Furthermore, these modes do not influence either the 
pitch or dumbbell motion; (c) because of the".’^rcsence of the dumbbell the natural 
frequencies and the mode shapes of the symmetric modes of the beam differ from 
those of the free-free beam. The frequencies and mode shapes can be easily ob- 
tained by replacing the dumbbell by a concentrated mass, equal to that of the 
dumbbell, at the center of the beara.^ 


^ STABILITY AJ^ALYSIS - BEAM WITH STABILIZING DUMBBELL" 

For small amplitude pitch and dumbbell oscillations and small deformations, Eqs. 
(9)-(ll) can be linearized to yield the following equations: 

e"+70'+(k-3)0-ca'-ka4-E(Ie'+ke » 0 . (12) 

n n n z 

a'’+c,,ca'+(c_k+3)a-c-c0’-c,k0-E(ce'+ke )c,C^”^ « 0 (13) 

i. I J. Xn nnxz 

e''+(n2-.3)E {k(o-e)+c(a'-0’)}C^^^J /M f2+Z('cE’+ke = 0. (14) 

n'n n ' 'zynm'mm'z 

where a j (m,n « 1,2 , . . .) and M « mass of the beam for all 

.4 z yzz'n*'’ ** n 


The small amplitude stability of the solutions to Eqs. (12) -(14) is governed by 
the roots of the characteristic equation. For the independent symmetric modes 
the system characteristic equation contains the separate factors 

Cs^+(n^-3)3 = 0; W =» 0) /m . (15) 

'n' *nn c 

for all the symmetric modes. For the remaining anti-symmetric modes the system 
characteristic equation can be expressed as / 


s^+cs+k-3 

-(cs+k ) 

(cs+k) 

• • (cs+k)C 

-(cs+k)c^ 

t 

s^+(cs+k) c j^+3 

• 

( 1 ) 

-(cs+k)c^C^ ^ 

• 

c » «> 

9 

(«4^)C^^*^/12 

# 

-(M+k)cf"Vl2 

In 

• • '^nn 

• 



• 


( 16 ) 


wiiere, \ 

■A » (cs+k)C^”®^ n^m 

A » s2+(7s+k)C^””^+ - 3 

nn ' ' z n 

If in a particular model only a finite number of elastic modes are considered in 
the model, then the associated characteristic equation can be easily obtained 
from Eq. (16) by deleting the rows and columns corresponding to the neglected modes. 
The stability analysis of a few selected truncated models is presented in the fol- 
lowing sections* . ' 

The Case of a ''Rigid Beam . 

Here it 1s assumed that = 0 for all n, and from Eq. (16) the characteristic 

equation can be obtained Dy using only those elements contained in the first two 
rows and the first two columns, as 

s** + (l+Cj^)cs3+(l+(:j^)ks2+3(l-c^)cs+3k(l-Cj^-9 - 0 (17) 

With an application of the Routh-Hurwitz criterion and noting that c. “ 
the following necessary and sufficient conditions for stability result: ^ 

c>0;k>0;c^<l;k> 3/(l-cp (18) 

* • 

The condition: c. < 1, implies that the dumbbell moment of inertia (X^) should 

be greater then tne-beam (pitch axis) moment of inertia. In the" limit as the 
spring stiffness (k) tends to infinity, the characteristic equation for the 
lower codes approaches: 

s2+3(l-C3_)/(l+cp - 0, (19) 

the characteristic equation the system would have if the dumbbell were to be 
rigidly connected to the beam at the center and pitching with it as a single 
rigid body. 


The Case of a Flexible Beam with only the First Anti-Symmetric 
Mode Included 

In this case it is assumed that 0 for all n 2 (n * 1, first elastic 

mode). The characteristic equation of this truncated model as obtained from 
Eq. (16) is given by: ' 

s^+cs+k-3 

-(cs+k) c^ 

(7s+k)c(2)/i2 

z 

where ^ 

After formal expansion and application of the Routh-Hurwitz necessary conditions 
for stability involving the sign of each of the coefficients in the characteristic 
equation, the following necessary conditions for stability must be satisfied: 


-(cs+k) 
s^+(cs+k) Cj^+3 
-(Is+k) 0^2^12 


(cs+k)C^^^ 
-(cs+k) c^C 

^,2 


( 2 ) 

z 


= 0 


( 20 ) 


+(cs+k)c(^>^^ + fl? 

Z JL 


3; J22 *= 
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c > 0; k > (3-n|)/a+Cj^+C^^»^^> 

^2 ^ ' k>9/{(l+Cj^)n|-6c^} 

> 3{C^^»^V(1-Cj^)+1) 

it > 9(Jl|-3)/{3(n|-3)a-Cj^)-9 (21) 


The Case of a Flexible Beam with More Than One but a Finite Number 
of Antl-Symmetric Modes In the Model 


In this case • 0, for n>m where m is a finite Integer greater than or equal 

to 4. As the number of modes retained in the model increases, the direct expansion 
of the detemd.nantal equation, (16), is algebraically complicated, and an alternate 
algorithm^ was used to determine the coefficients of the characteristic equation. 
For the implementation of this algorithm, it is necessary to rewrite Eqs. (12) - 
(14) in the standard state variable form 


X' - AX 




( 22 ) 


f 

where X’“ [0,a, e,,... e ,. 0',a', e, * , ... e* ]' 
1> m ^ 1 m 


The characteristic equation obtained by this algorithm can be solved for the 
characteristic roots by any of the many algorithms, for determining the roots of an 
algebraic polynomial equation. (In this case an algorithm based on the modified 
Bairstow®»^ method was used.) 


NUMERICAL RESULTS - BEAM WITH STABILIZING DUMBBELL 

Fijgure 2 depicts the root loci of the least damped mode with the spring stiffness 
(k) as the varying parameter for increasingly complex models of the beam-dumbbell 
system. The root loci are symmetric about the real axis and only the portion above 
the real axis is shown. The upper curve labeled corresponds to the case of 

the rigid beam. It can be noted that as the spring’^stiffness increases the_roots 
tend toward the imaginary axis. For the rigid beam case, from Ineq. (18), k‘ >30 
to ensure stability for the value of c^= selected. 

When the first anti-symmetric mode is Included in the model (solid curves in 
Fig. 2) the locus is obtained by assuming a particular value for Tlie 

following observations can be made: 

(a) with increasing values of the spring stiffness , k, the roots of the 
characteristic equation move toward the imaginary axis^ 

(b) the addition of an elastic mode into the model has the effect of 
moving the eatire root locus toward the origin; 

(c) the minimum value of spring stiffness required for stability is 
increased. 

With the further inclusion of the first two and the first three anti-symmetric 
modes of the beam into the model, it can be noted from Fig. 2 that a further 
slight deterioration of the stability characteristics results. The minimum 
value of hinge spring stiffness required for stability is further increased in 
these cases. 
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A representative transient response to an initial displacement about the local 
horizontal is »?hovm in Figure 3. Tlie long time required to completely damp 
the system oscillations in the lowcast frequency mode may be attributed to the 
fact that the dumbbell inertia is only slichtly larger than the beam's pitch 
axis inertia in this example (cj^«0.9). Increasing the dumbbell inertia can be 
accomplished by Increasing the connecting linkage length within certain limits 
and/or increasing the tip masses at the expense of the overall useful payload. 

It is thought that the use of an active control system together with the hinged 
spring damper may prove to be the most efficient way to remove transient motions 
and also to achieve shape control in all modes. 


A THIN, UNIFORM FLAT PLATE IN ORBIT WITH ITS MAJOR SURFACE 
NOMINALLY IN THE LOCAL HORIZONTAL PLANE 


It is’ known that a flat plate with its surface normal nominally along the local 
vei’tlcal is gravitationally unstable in the absence of external restoring torques. 
In order to overcome the destabilizing effect of the gravity-gradient torque on 
the plate, one can either apply active control torques or adjust the moment of 
inertia distribution of the system such that the gravity-gradient torque now 
becomes stabilizing. As in the case of the^iseam, the plate also can be gravi- 
tationally stabilized by attaching a rigid, light weight dumbbell of proper 
moment of inertia at the center of mass of the plate (Fig. 4). The dumbbell is 
assumed to be attached to the plate by a spring loaded double gimballed joint. 

Thus the dumbbell possess two degrees of freedom described by the angles y and 5 
as shown in Fig. 4. Damping may be assumed to be present at the hinges of the 
gimball. 

The following assumptions are made in deriving the equations of motion for this 
system shown in Fig. 4: (a) the plate has a constant thickness which is much 

less than the other characteristic dimensions of the plate. The mass and stiffness 
properties are uniformly distributed throughout the plate; (b) the center of 
mass of the plate is moving along a circular orbit in a spherically symmetric 
gravitational field of the earth; (c) the elastic displacements in the plane 
of the plate are negligible compared to those normal to the plane of the plate; 

(d) the plate is completely free (unconstrained) ; (e) the mass of the rigid 
link connecting the two tip masses of the dumbbell is negligble compared to the 
tip masses; (f) the center of mass of the dumbbell coincides , with the center 
of mass of the undeformed plate; (g) the attitude angles and vibration ampli- 
tudes are small; (h) there are no external disturbance and control torques. 


With the above assumptions one can derive the linearized equations of motion 
in the following nondimensional form: 




9" - 36-c y'-k y+E(c e’+k e )C 
y' y* n y n y n z 


(n) 


= 0 


e" + («2-3)e -{(c y»+k y)J C^’^^+(c 6'+k 6)J 
n ' n n y y ' y z z z z y n 


,(n) 


(23) 

(24) 

(25) 


+E{(c e’+k e 
m y m y m z 


+(c e'+k e )C 
- ~ z m 


t. 

z 'm 


(mn) 

y 


) = 0 


(26) 


/ 
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<a?) 


Y"+’cy(l+c^M ** ® 

6'*+ 7^a+C2)6’+{4+kj(l+C2)}^+8«-2iji’-(l+C2)g(CjjeH^^^ * 0 <28) 

where, i|>, «J>, 0 ■ yaw, roll and pitch angles, respectively ; ^ "* 

amplitude; M ■• n**^ modal mass; Z * characteristic length (eg. radius of a cir- 
cular plate or length of the side of a square plate); Y» <5 * dumbbell angles; 

»* principal moments on inertia of the plate; 1^ "• principal moment of 
inertia of tha dumbbell; Qjj"* /^•x » “ torsional spring stiffness ; 

Cy,C 7 > ■■ damping co-efficients; and, 


1- 

V*d» =2 ■ V^d' \ ■ 


• ]t M 

' z 

k /J 
z' z c 

m 

c,./ J . 0 ) ; c ■ c / J u ; 

0<"> . 2 

gx(n) 

1 

y 

y y c* z z' z c' 

y 

sy 

'y ■“ 0, z ' 

•a t 



cf™>- , 

J 


z 

ly - 0, z - 0 

y 

y y 

y ' m 

,(mn) 

'z 

- J 

z z z m 



t * 


A study of Eqs. (23) through (28) reveal that: (i) in general, the pitch, roll, 

yaw, dumbbell angles (y and 6) and elastic motions of the plate, are all coupled 
to each other; (ii) the elastic modes for which C^, - “ 0 and )=» 0 can neither 

influence \|»,<{),0,y ^ aor be influenced by them? however, their frequencies , and 

mode shapes are modified by the dumbbell; (ill) the elastic modes for which 0 

and 0 are not influenced by and 6. However, they are directly coupled to 

0 and’^Y; (iv). the elastic modes for which 0 and 0 are not influenced by 

0 and y; however, they are directly coupled^ to and 6; (v) the elastic modes for 

which neither nor vanish couple ij>,<{>,0,Y and 6 motions; (vi) the natural 

frequencies and^raode shapes of the elastic modes for which the hinge point (y»0,z*0) 
lies on a nodal line, are not affected by the dumbbell. 

The stability of the solutions of the Eqs. (23) -(28) to small initial conditions 
depend on the characteristic roots of the system. Theoretically, the characteristic 
determinant of Eqs. (23) -(28) is of infinite order. However,, if only a finite num- 
ber of elastic modes of the plate are retained in the model, the characteristic 
determinant corresponding to the truncated model is of finite order. A stability 
analysis of a few such truncated models is discussed in the following section. 


STABILITY ANALYSIS-PLATE WITH STABILIZING DUMBBELL 
The Case of a RiRld Plate 

For a rigid plate, 0.. for all n.‘ Hence, 0 and 0. This results in 

(i'»,(|),6) and (9 ,y) being independent of each Other. The characteristic roots of the 
system Eqs. (23) -(28) for this case are given by the roots of the following equa- 
tions : 
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Roll-yaw and 5 motion 


8^-n 


0 

2s 

(s^-4) 

-(c^s+kp 

-2s 

8 



™ 0 (29) 


Pitch - y rootion 


a‘*-tay(l+c^)3 3+ky(l+Cj^)s^+3Cy(X-Cj^)8+3ky(l-Cj^)-9 - 0 

where, ■« s^+Cl+c^) (Cj^s+kjj)+4 (30) 

From Eq. (29) it can be observed that for plates with, 0^ "O (e.g. square plates, 
circular plates etc.), the roll-yaw raption tends to be unstable, which is indi- 
cated by the double root at the origin of the complex plane. This instability 
in roll-yaw may be overcome by introducing a slight asymmetry into the plate 
such that, |njjl«l but not equal to zero. This may be achieved by placing 
small concentrated masses at suitable locations on the plate as dictated by 
the following necessary conditions for stability which are obtained from Eqs. 
(29)-(30). 


c >0 ; c^>0 ; Cj^<l ; C2<(4-J2^)/(2-Q,^) ; fi^<0 ; k >3/(l-Cj^) ; 


Hence, the slight asymmetry Introduced by the praall *eoncentrated masses should 
be such that, n„<0, i.e. . . 

.-JC .-tf Z y ^ 

It is Interesting to note that in the limit as k -»• « and k «, the charac- 
teristic roots corresponding to the low frequency modes of the system tend 
toward the roots of the following equations, 


s‘^+{(l+c^)(2-fi^)+2}s2/(l+c2)+4ftj^(c2-l)/(l+C2) “ 0 ' (31) 

and s^+3(l-Cj^)/(l+Cj^) 0 (32) 

Fig. 5 shows f.he locus of the roots corresponding to the least damped_mode of 
Eq. (29) with Ic^ as the parameter along each curve. As the value of k^ is 
Increased, the roots move toward the imaginary axis to the value as given by 
Eq. (31). Also, we note that with higher values of |Hx 1» roll-yaw sta- 
bility, of the system is improved. The locus of the roots corresponding to 
the least damped mode_of Eq. (30) is shown by the solid curve in Fig. 6. Here 
also we note that as ky is Increased, the roots move toward the imaginary axis 
to the value as given by Eq, (32) . The effects of flexibility on the system 
stability characteristics is considered in the next section. 


The Case of a Flexible Circular Plate 


In order to investigate the effects of flexibility on the dynamics , we have 
considered a circular plate in the present paper. It is well known that the 
elastic mode shapes of a thin, uniform completely free circular plate are 
characterized by a nodal pattern which consists of nodal diameters and con- 
centric nodal circles centered at the center of the plate. Mathematically 
these mode shapes are given by,®»^ 


.(n) 


X 
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(33) 


where » p 

j 




D - 


h w 


5 - 


J ,I 
P P 

A ■ • G. 

j.p j»p 


Based on Eq< 


number of nodal diameCers (p » 1,2,...) 

number of concentric nodal circles (j *• 1,2,. ».) 

the frequency parameter, (w. ; 

J »P 

natural frequency 
mass density 

flexural rigidity, Eh3/12(l-v^) 
thickness of the plate; v * Poisson’s ratio 
Young’s modulus 

nondimenslonalized radial distance from the center rf 
the plate (0j<55 1) 

It ‘ , _ 

polar angle measured with respect to 'the y axis 
arbitrary phase angle 

Bessel and modified Bessel functions of the first kind 

constants whose value can be adjusted to yield the modal mass 
equal to the mass of the plate® 

(33) one can easily show that for circular plates, 


\^0 for P5i*l 



for p - 1 

\^0, for P^l 


(3A) 


(35) 


Thus, Eqs. (23) through (28) are coupled to each other through the clastic 
modes which have only one nodal diameter. The other elastic modes are inde- 
pendent of 


In the previous section it was assumed that small concentrated masses were placed 
on the plate such that, to prevent" the instability in the roll-yaw degree of 

freedom. The concentrated masses arc assumed to be placed at the ends of the 
diameter along the roll axis (z-axis) of the plate. Because of this asymmetry 
in the plate, the arbitrary phase angle, 3o> i>^ Eq. (33) can now have the values 
0 or ^/2 only A® This implies that the nodal diameters can have only one of 
the following two orientations: (1) a nodal diameter along the roll or z axis 

of the plate (3o®0) » and (2) a nodal diameter along the pitch or y-axis of 
the plate (0 q=it'/ 2). The natural frequency and the mode shapes of the plate are 
not affected by the concentrated masses for the first orientation of the nodal 
diameter. However, the small concentrated masses slightly perturb the values 
of the natural frequencies and the mode shapes of those elastic inodes for which 
the position of the concentrated masses do not lie on a nodal line. 
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To a first order appro 2 ?iraation , the values of the perturbed natural frequency 
and the perturbed mode shapes are given by,^® 


P n n 


’’ ■ (ni*.!) 


(37) 


where, u - perturbed natural frequency and mode shape 


P P 




X 

6M - 


xmperturbed! natural frequency and rocde shape •. 


total mass added 


« -a. 


m 


m2 .,m2 

m n 


n • 




In the following sections the dynamics of a flexlbe circular plate for the 
previously mentioned two possible orientations of the nodal diameter are 
discussed separately. 

• 

The case of a nodal diameter along the pitch axis or y axis of the plate (3n°IT/2) , • 
For this orientation of the nodal diameter for all n (E'q. (34)‘. Hence, 

(^,(}>,d) motions become completely independent of pitch, y and elastic motions of 
the plate and thus the plate behaves like a rigid body in the roll-yaw degree of 
freedom. However, the 0 and y motions are coupled to the elastic motions of the 
plate through the terms containing c(’^) in Eqs. (25)-(27). The stability of the 
so^.i^tlons of these equations to small initial conditions is governed by the roots 
of i characteristic equation which theoretically is an infinite degree polyno- 
mial equation. However, in practice since only a finite number of modes are 
retained in the mathematical model, the degree of this polynomial equation will 
be finite. A stability analysis of such truncated models is presented in Fig. 6. 

It can be noted from Fig. 6 that the characteristic roots corresponding to the 
least damped mode, of the system behave similar to those discussed in Fig. 3. 

Tlie case of a nodal diameter along the roll axis or z axis of the plate (g^«0). 

For this orientation of the nodal diameter 0. for. all n (Eq. (35)). Hence, 

(0,y) motions become completely independent of (J>, 5 and elastic motions of 

the plate and thus the plate behaves like a rigid body in the pitch degree of 
freedom. However, i|i, 6 are coupled to the elastic motions of the plate 
through the terms containing in Eqs. (23), (24), (26) and (28) v A para- 

metric study of the least damped mode characteristics for these equations re- , 
suits in a root locus plot similar to rhose presented in Figs. 2 and 6. 


CONCLUSIONS 

Space ’structures with their maximum moment of inertia axis along the local ver- 
tical are gravitationally unstable. Such structures may be gravitationally sta- 
bilized by attaching a ligtit weight dumbbell with heavy tip masses. 


3-12 


An analysis of the dynamics of two such systems, a thin fle:dble beam nominally 
along the local horizontal and a flat plate with its surface normal nominally 
along the local vertical, is presented in this paper. 

In the case of the flexible beam, the dumbbell, motion excites only the anti-> 
symmetric elastic modes of, the- beam.- .However, the frequencies and the mode 
shapes of the symmetric modes of the beam are modified by the presence of the 
dumbbell. 

For case of flat circular plates, the dumbbell can excite only those elastic 
modes with only one nodal diameter. Also, the dumbbell modifies the value of 
natural frequency and mode shapes of axi-symmetric modes ci the plate. The 
small concentrated masses added to the plate for roll-yaw stability results 
in only two possible orientations for the nodal diameters. These concentrated 
masses slightly perturb the value of the natural frequency and mode shapes, 
if they do not lie on a nodal line . . 

For all the cases presented in this paper, the characteristic roots corresponding 
to ^he least damped mode of the system move toward the imaginary axis as the 
torsional spring stiffness at the hinge is Increased. The stability of the 
truncated model deteriorates with the increase in the number of elastic modes 
retained in the model. 

In order to damp the motion of the system in all its modes , especially the low 
frequency modes, the use of active dampers (control systems) are needed. . However, 
it is thought that with the use of the passive gimballed dumbbell stabilization 
device together with active controllers, both the peak forces ns well as fuel 
consumption could be -significantly reduced. Such a. study would 'represent a logical 
extension of the pres;:jnt paper.. 
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Fig. I:’ • Dumbbell Stabilized Flexible Beam With Nominal 
Orientation Along' Local Horizontal 
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Fig. 3: Time Response of Dumbbell Steblllsed Flexible Bean la Orbit 
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Fig. 4': Flexible Circular Plate Nominally in the Local 

-Horizontal Plane with the Stabilizing Dumbbell. * 
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Fig. 5: Root Locus of Least Damped Mode —Rigid Circular Plate 

with Dumbbell 
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Chapter 4 


ON THE MOTION OF A FLEXIBLE SHALLOW SPHERICAL 
SHELL IN ORBIT 

In this chapter an analysis of the motion of a thin, shallow spherical 
shell type structure in orbit nominally pointing towards the earth 
along the local vertical is presented. It is known that in this 
orientation, the system is gravitationally unstable. Such structures 
may be passively stabilized with a connecting dumbbell. Hence an 
analysis of the motion of a shell structure with a stabilizing dumb- 
bell is also presented in this chapter. 


DEVELOPMENT OF EQUATIONS OF MOTION-SHELL AXIS 
ORIENTED ALONG THE LOCAL VERTICAL 

Fig. 4.1a shows a shallow spherical shell with the various notations 
defined. The assumptions that vete made in deriving the model for the 
shell are as follows ; 

(i) the mass and elastic properties are distributed continuously 
and uniformly throughout the domain of the shell; 

(ii) the thickness of the shell is small as compared to the height 
of the shell; 

(iii) the ratio, height (H)/base radius (A) is much less than 
unity (condition for shallowness) ; 

(iv) the edge of the shell is completely free; 

(v) the elastic deformations perpendicular to the symmetry axis 
(i.e. x-axis) of the shell are negligible compared to the de- 
foxnnations parallel to the symmetry axis, i.e. only transverse 
vibrations are considered; 

(vi) the center of mass of the shell is moving in a circular orbit. 


Tlie transverse vibrational mode shapes of shallow shells can be conveniently 
expressed in a cylindrical system of co-ordinates (r ,3,x ) can be related 
to the X ,y ,z system (see Fig. 4.1a).^'^ ^ 


In order to evaluate the coupling tierm^in Eqs. (2.5) and (2.14) for the 
present case, let us express, r^ = ^2 ('^•l) 
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where , 


r, ■ the vector from th;i center of mass of the shell 

^ to the origin of (x y ,z ) 

0 0 0 

“ the vector from the origin of (x y z ) to a generic 
point on the shell. ^ 


Since the shell is assumed to be-.Qompletely free, by. Eqs. (2.16) and (2.17) 
we have, 

/ q dm « 0 

V 

/ q dm =“ 0 (4.2) 

V 

00 

f rQXq dm 0 


Using 


n 

•Si, 

zcp 

ra Jmn 


Eqs. 

Q(n) 


(4.1) and (4.2) in Eqs. (2.5)- (2.14) one can easily show that. 


» / [2r2x(u)xq) + r2x(o)xq) + qx(Wxr2) ~ (r2*w)(a)xq) 


-(q*w) (wxr 2 )] dm 

(4.3) 

= / (r 2 xM^‘^^q+qxM^‘^^r‘ 2 ) dm 

(4.4) 

v 

• 


a / • 0 )xrA+<j"^^^ •tox(o)xr«)] dm 

(4.5) 

a f •O)x(ti)xq) dm 

(4.6) 

V 


» / dm 

V 

(4.7) 


Sg 

m ran 


= / dm 


f(c) 


T-MT,"’ 
4 4 


where, M 

=• transformation matrix given by Eq. (A. 2) 
M = matrix operator given by Eq. (A. 4) 


(4.8) 


In deriving Eq. (4.6) use has been made of the assumption that th^e^elastic 
deformations are parallel to the symmetry axis only, i.e. q = Iqli* 
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The following observations can be made from Eqs, (4.4)-(4.8); 

(a) the rigid body rotational mode and elastic modes of the 
shell are coupled to each other by both inertia and gravity; 

(b) the earth's gravitational field can excite elastic motion in 
the shell through rigid body coupling terms. 


After expressing the various vectors in Eqs. (4. 4) -(4. 8) in the 

cylindrical co-ordinate system (r ) we have, 

o c 


r 


2 


A A 

r e +x i 
c r c 


(4.9) 


— " '' 
tu ■ 0) e -HDoeo-KO i 
r r 8 3 X 


AAA 


W^i-H0yj+0)^k 

(4.10) 


(4.11) 

X 

(4.12) 


By substituting Eqs. (4.9)-(4.12) into Eqs. (4.3)-(4.8) one obtains 

q(o) ^ y (dJ-O) i) - r 0) i} 

V n ^ c' X ' c r 

i ^ , A 

+ A {x (oj-o) i) - r 0) 1} 
n c X ' c r 

-L • ^ 

+ A {x (oJ-O) i) - r 0) e } 
n c' X c X r 

A A 

- A (r 0 ) +x 03 )(0Jj-O3 k) 
n c r c X zr y ' 



- A 0) {o)qX e -(03 X -r 03 )eQ-r Woi}] 
nx pcr rc cx 3 c3 '^^x 

(4.13 

5<n) 





(4.14) 


=» - / {r (Lq-Hx (0)^-Ko^)-r 03^03 } t()^’^^dm 
c y z' cxr ^x 

(4.15 


= -A (o3^-ko^)S M 
my z mn n 

(4.16) 

Sn 


(4.17) 
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(4.18) 




where, are elements of the synunetric matrix, given by, 

’ ij 


_ „ . Symmetric 


M22C^3-M32s23+M2^s^0 I 

M2j^s3+M3j^c8 |(M22-M23>8234M32c23 M228^3+M32^23 

^ +Mj3=^6 _ 

(c' 

The other elements of the matrix are generated by noting, 
are the elements of the matrix M (ref. 3) . 

After substitution of the following 

Sj. “ c3j + s3k 

e^ * -’s3j + c6k 

03 » W c3 + 0 ) s3 
r y z 

03o ■ -OJ s3 + 0 ) c3 
3 y z 

into Eqs. (4.13) - (4.18), there results 

+{2(A„yA„J,^)lf ) + 

+(03W)I^''^} j 
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+ {2M21I4" - 2M3j^l5 -M23I7 +2 ^ 

(4.22) 

<<? 

’ n 

■ Vs"’- 

(4.23) 

'ran 

- -\<“>z>«,nn “n 

(4.24) 

8n 

^I’^l ^^1^2 ^ 3 r 3 

(4.25) 

^ran . 

\^ll'^mn*n 

(4.26) 

where , 

/ r^c3tl)^“^dm; / r^sg(|)^''^dm; 

! X -.2p(t)^'^^dm; / x s23(|)^^^dm; / r c30 dm; 



/ ^£.83$ dm 


Since, the shell is assumed to be completely free, the mode shapes satisfy 
the condition in Eq. (2.17). 

i.e. 

; r„x$^''^dm » 0 

V ^ 


With r^ 

» 7 +r and Eq. (2.17) translates into, 

12 X 


» 

/ r„x(j)^’^^i dm * 0 
2 ^x 

V 


d. • 6 • 

- » 0 


! 

1 
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Hence » 


(4.27) 


i<")h 0 


(4.28) 


The mode shapes of a completely free shallow spherical shell are charac- 
terized by a nodal pattern consisting of a set of j concentric nodal circles 
centered on the symmetry axis and p nodal diameters. The mode shape as 
given by Eq. (B.4) in Appendix - B is in the form 


cosp(a+6^) 


(4.29) 


where, ? ■ r /A and I ■ the base radius of the shell . 
o 

Using the. form of in Eq. (4.29) we can evaluate the integrals, 

through ly ^ as follows: 

- / X ^^”^dm 

1 V c * 


p X (?)W. (K) cosp(8+fi ) ?d?d3 

0 0 0 


i • G • I 


2m 6 

1 op X. 


(4.30) 


where, /^x^(?)W, 


*0 0 

m ■ mass of the shell, irpH^ 


/ 1» P ■ 0 
^ 0, p ,4 0 


Now consider, 


/ X c 23 (l>^'^^dm 

H y C X 

" /^/^'^x (?)W. (?)cosp(8+3rt)cos23 ?di;d3 

0 0 J’P 


m cospBq 


Similarly, 


-m stapSp «2p4f 


(4.31) 


(4.32) 
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(4.33) 


I<">- mi cospPo 

- mi slnpBg 


«here, j(j). />* ( 5 )» <;)Cd5 

^2 0 ° 


l(j)« (5) dC 

0 

V ^oipjk <“ 


(4.34) 


The integrals 1^2^ and are evaluated in Appendix C. Their values 

MWA A<l «*A«% W«« ^ 


are given by, 


(4.35) 


" OP 9 1-1 ovniT^ o {2X, ,(I,(X, J+1)-8I,(X, J} 


X2 2R ‘j,2 ^12RDXJ^2 J»2 ,2"*' ''‘l'“j,2' 




(4.36) 


“ ^rr: "j.a ^«-«"0<^3,3>-^S,3"l<N,3>-^d,3^2<S.3» 


j,3 


D. 


■'■ i\ , ^^j.3^2^^j,3^”‘^^j,3^1^^j,3^’^®^0^^j,3^"®^ 

j » -> 


(4.37) 


where, V = Poisson’s ratio 


A ,C ,D = constants in the mode shape function in Eq. (C.4) 
3,P d,P 1,P , 1,2,...) 

After substituting Eqs. (4.21)-(4.28) into Eqs. (2.2) and (2.3), we 
arrive at Eqs. (4.38) - (4.41) for the motion of the shallow shell 
£n orbit. 
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Vaw: u) +Q W u) 1L„+G„/J 

X X y z X 73 x' X 


Pitch: w -HI 0 ) 0 ) +2t(k w '»*A w -A - 

\f \r V n ' r% If ii \f v v ' T 


x<"> 


(4.38) 


y y z X n’ n y n y u z x' 


y i 

I<«) 

Roll: u)^+n^u^y2E(A„u^+A_^io^+Au^Uy)^ 

■ z 

"^z^l2'*'n^ 2^^’22’"S3^^6 ^'*'^^^2l''4 ^““^^^’31^5 

# JS 

Modes: A 4«^A '(<oW)A - M^ -jf- + M A„+ jf- 

n n n 

where, S7^ - (J^-Jy)/J 3 , : Sly - ' % “ <'’y''^x^ > 

“w/w ;e «A/il;Jl“ base j^adius of the shell; 

In order to examine the stability of the system response to small initial 
conditions, Eqs. (4.38)“(4.41) are linearized assuming- small amplitude 
pitch (0) , roll (<i) , yaw (i()) and clastic displacements (A^) and denoting 

|;^ ( ) « ( )V; t « w^t. 

As a result of this we arrive at the following linear equations of 
motion for the shallow shell; 


(n) 


(n). 


. w A C 
(rt) N \ ^ ^ JL 



(4.42) 

(})"+4n (|)+(1-S^ )i|;’ » C /J 

Z 2 Z 2 C 

(4.43) 

0"-3Q 0-2Ee'i^”^Vj “C /J 

y n n 1 ' y y y c 

(4.44) 

e"+(a^-3)e +20 ’I^"Vm /M 

n n n I n 1 ' n n n c 

(4.45) 


An examination of Eqs. (4.42) - (4.45) reveals the following points 
(external disturbance assumed to be zero): 

• in the linear range of operation, the motion in the roll-yaw 
degree of freedom can be studied independently of the pitch and elastic 
motions; 

• the pitch and elastic motions are coupled directly to each other 
through their races; 

• s;l ce 0 for all elastic modes except for the axi-symmetric 

modes (i.c. modes with no nodal diameters), only axi-symmetric modes 
are responsible for the coupling of pitch and the elastic motions, Non- 
axisymmetric modes are independent of the rigid body rotational motions, 

wand 6; 


(4.39) 

(4.40) 

(4.41) 


4-8 


• the axieywnietric elastic modes are subjected to a constant 
excitation force due to the orbital motion and gravity effects; 

• the pitch and the roll-yaw motions are unstable about the 
present nominal orientation because of the unfavorable moment of the 
inertia dlstri.^ntion. 

In order to stabilize the pitch motion, a passive stabilization 
procedure using a light weight dumbbell, as proposed in the case of 
flat plate, is also considered in the next section. 


GRAVITATIONALLY STABILIZED SHALLOW SHELL IN ORBIT 


A gravitationally stabilized shallow spherical shell in a circular orbit 
is shown in Fig. 4.1b. Tlie stabilizing dumbbell is assumed to be hinged 
to the shell at its apex by a spring loaded double gimbal joint. Thus, 
the dumbbell has two degrees of freedom with respect to the shell. Damp- 
ing is assumed at the hinges of the gimbal. • The two tip masses of the 
dumbbell are assumed to be connected by a rigid link of negligible mass. 


Tlie reaction torques on the shell due to the relative motion between 
the dumbbell and the shell are given by, 


{k^(6-a^)-h:^(6-a^)> ay 

(4.46) 

ky(Y-Oy)+Cy(Y-<?y) 

(4.47) 

{k^(5-a^)+c^(5-0^)} cY 

(4.48) 


where, k ,k * torsional spring stiffness at the hinges of 
^ ^ the gimbal 

c ,c = damping co-efficients at the hinges 

y 2 

Y,6 =* dumbbell angle defined in Fig. 4.1b 


0 ,cr 
y z 


small rotations of the shell at the gimbai 
due to elastic deformations. 


The modal components of the reaction torques in Eqs. (4 . 46) -(4. 48) can 
be expressed as, 


E 

n 


G 


y 


3z 


y=0 

z-0 


+ C 

z 


9y 


y=0 

z=0 


(4.49) 
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Eqs. (4.38)-(4.41) together with eqs. (4.46) ^(4.49) describe the motion 
of! the shallow shell. The equation oi! motion of the dumbbell in the 
TjCxdyd^d) frame (see Appendix A) are given by, 




d "d ‘'d 


“z/'V.V ■ 

d ■^d d d 


(4.50) 

(4.51) 


where, U) ,U) ,co ■ angular velocity components of the dumbbell, 
’^d ^d ^d 

G ,G “ gravity-gradient torques on the dumbbell. 

^d ^d 


In order to Investigate the stability of the system about its nominal 
local horizontal orientation, Eqs. (4.38)-(4.41) , (4.50) and (4.51) are 
linearized for small pitch (0), roll (({)), yaw (ijj) , Y» elastic dis- 

placements. Tlae linearized equations are presented in Eqs. (4.54)-(4.59) . 
In deriving these equations we have used the results , 


a 

y 


8z 


y=»0 

z“«0 


E E 

n n z 


a 

z 


9q 

^x 

“ 

3y 


y»0 

z^O 


E e 

n n y 


where. 


e =■ A/il ; 
n n 


C<n) 

y 


a 


9(|), 


(n) 


X 


9y 


y=0 

z=0 



9z 


y=0 

z=0 


(4.52) 


(4.53) 


Since, C = 0 for all modes except for the modes with only one 

nodal diXraeter^ it is evident from Eqs, (4. 54) -(4. 59) that in the linear 
range the pitch* roll, yaw, y and 6 motions are coupled .to the elastic 
motion only through the axi-symmetric modes for which iivf 0, and the 
modes with only one nodal diameter for which 0, Q. 


r - )‘i>' = 0 

A A 

(j)" + 4S2 (j)+(l-Q )i|i’ = c 6'+k 6-E(c e’+k e ) 

^ z z ^ z z n r, n z n y 


(4.54) 

(4.55) 
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(4.56) 


0" - 0 - - a Y' -I- k Y - S(c e'+k e )G 


(tv) 


m S»» Aim w. 

n n 1 J 

j,(n> 

en H- (n^-3)e_,+20' 5^. 

n 


y 

(n) 


11 y n y n z 


HT- + <=yY-+u/> mV ^ 

11 n 

J 


+ (c 6’+k 6) - 2 (c E'+k e 


m ' y n y n z 


,(inn) 


- 2(c e'+k E )C™(n - 1,2,...) 
in z n z n y 


Y" + Cy(l+cpY'+ {3+ky(l+c^)} Y + 3(My)0 - 

+ 2 2e’ - 0 

n n 1 I„ 


(4.57) 


(4.58) 


5" 4-7 (l+c-)6' + {4+k (l+c,)}6 + 4(1-J2 )i> - (1-Q )\p' 
z ^ z ^ z z 


- (l+c„)2(c £’+k E « 0 
z n z n z n y 


(4.59) 


In the following sections the stability of the system when a finite number 
of ‘elastic modes are retained in the model is considered. 

. < t \ • . 

The Case of a Rigid Shallow Spherical Shell 

For rigid shallow shells , 0 for all n. This results in decoupling 

of Eqs. (4.54), (4.55) and*(4.59) from Eqs. (4.56) and (4.58). Thus, the 
stability of the system for small initial conditions is governed by the 
following two independent characteristic equations. 


S^~n -'■vl+fi )s 0 

X x' 

(1-fi )s s^4-4fi “(c s+k ) 

z z z z 

-d-J^ps 4(14^^) 


= 0 


(4.60) 


s^-3fl -(c s4-k ) 

y y y 

3(l-K^y) 

where, A^^ = s^4-(l4-c,) (c^s4-k^) 4- 4 
A^^ = s^4-(14^^) (c^s4-k^) + 3 


(4.61) 
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For a perfectly symmetrical spherical shell, - 0, This results in 
Eq. (4.60) having a double root at the origin and, thus, rendering the 
system unstable in roll and yaw. To prevent this instability, it may 
be desirable to introduce a slight asymmetry in the structure such that 
0<|J^^|«1. 

After formal expansion of the characteristic determinants in Eqa. (4.60) 
and (4.61), we arrive at Eqs. (4.62) and (4.63). 

■ 0 (4.62) 

s'*4ct^s^d<tgS^+agS4a^Q " 0 (4.63) 

where, “ c^(l+C 2 ) 

“ 5 + 

»3 ’ + 4-n^} 

+ 4{l+3n -2Q f2 } 

Z z X 

tte * -4c (c«n +1)Q 

^ Z Z Z X 

“6 “ 

“7 ■ 

a„ - r a+o,)+ 3(l-n ) 

0 y 1 y 

ci^0= 3ky (1-c^Qy)- 

-1 < Q < 0 and 0 < J2 < 1 
- z y - 

After applying the Routh-Hurwitz criterion for stability, we arrive at the 
following necessary and sufficient conditions for the stability: 

c >0 
z 

K ^ (n (^3)-5}/(1 +c^) 

Z Z X z 

Co <(4-5^^)/(f^o^x"^^z"^^ 

k >-4{l-f3S2 -2n }/{c„(l+3n^-fiX.)+ 

Z ZZXZ ZZX X 



n < 0; k > vIl 

X ’ z l+c„Si_ 

7 > 0; k > -3(l-J2J/(l+c.) 

y / / **• 

c, < l/fi ; k > 3Q /(1-c.n ) 

1 y y y 1 y 

and > 0 (i » 1, ..., 5) 

where, A^j^’s are the principal minors of the determinant 


“2 “l 


0 0 


“5 “4 “3 *'2 “l 

0 cij Cj 

0 0 0 a, Oj 


Tlie condition 0 implies that Jy 
two small concentrated masses at the 


y > J . This can be achieved by adding 
he enis of the diameter along the z-axis. 


As k and k become very large, the characteristic roots corresponding 
to tKe lower frequency modes are, to a very good approximation, given 
by the roots of the following algebraic equations 

s^-kx. ■ 0 ' 

4 6 

s^+cx “ 0 ^ 

s tu^Q u 

where, a, =* {c„(l+3f2 -f2 ^ ) 4- 4 - }/(l+c„) 

q Z Z Z X X z 

ttg ® -4(c2J^2,+l)fi^/(l+C2) 

3(l-Cj^ny)/(l+c^) 

The roots of Eqs. (4.64) and (4.65) are given by, 

^1.2 = ± 

"3.4 


®5,6 = ± “lO 
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0.9, 


For the values oi the parameters given in Fig. 4.2 i.e. c„ » 

* -.0001, * 0.9, H/a " 0.1 we have 

Si 2 “ ± i 4.2640143x10“^ 

s., , - + i 1.076055 
3,4 — 

3- . - + i 0.4090303 
5,6 — 

Fig. 4.2 shows the locus of the roots corresponding to the lowest frequency 
mode in the roll-yaw degree of freedom. The spring stiffness k is varied 
along each of these curves. It can be seen from these curves tSat by 
placing heavier concentrated masses at the ends of the diameter the roll- 
yaw stability can be improved. 


The Casa of a Shallow Spherical Shell with a Finite Number 
of Elastic Modes in the Model. 


Tlie mode shapes of a shallow spherical shell are characterized by a nodal 
pattern consisting of p (p =» 1,2,...) nodal diameters and j (j =* 1,2,...) 
concentric nodal circles centered about the axis of symmetry. For a per- . 
fectly symmetrical shell the position of the nodal diameter is arbitrary 
which is indicated by the arbitrariness of $ in Eq. (B.4) in Appendix - B. 
However, the addition of two small concentrated masses on the z-axis to 
make the shell slightly asymmetrical for roll-yaw stability, permits only 
two possible orientations of the nodal diameter: 

a. nodal diameter aligned along the roll axis (z axis) 

i.e. " 0. 

b. nodal diameter aligned along the pitch axis (y axis) of the 
plate i.e. 3 q * tt/2. 

Hence, the uncontrolled dynamics of the shell for each of these possible 
orientations of the nodal diameter is presented separately in the following 
sections. 


A. The case of the nodal didmeter along the. pitch axis (y axis). For 
this position of the nodal diameter, 3 q “ 'n'/2, in Eq. (B.4) in Appendix B. 

Hence, > s 

e''"'’ = 0 for all n 

y 


and 


cCn) 

z 


I 0, 


,1 ^j,l 

p 1 


P“1 


( 4 . 66 ) 


Hence, the roll (({)), yaw (ip) and 5 motions completely decouple from pitch, 
Y and elastic motions, i.e. (]jj, (j), 6) and (0, y, ......... E^) are two 

independent sets. Further, we consider only axi-symmetric modes and modes 
with only one nodal diameter in the present analysis. The other elastic 
modes are all independent in the uncontrolled dynamics . 
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The characteristic equation corresponding to the set (P, 6) is given 
be Eq. (4.60). Tiujs, the shell behaves like a rigid body in the ip, <p, 

6 degrees of freedom. The characteristic equation corresponding to the 
set (9, Y, is given by, 


s^-3ny 


3(14fly) 

V 

2I^^^s/M^A 

• 

1 

• 

2I.^“^s/M Z 
1 n 

0 

n 


-d+Cj^) (l/Jy)'. . 



• • • 


(4.67) 


where, __ _ 

=■ s^+(l+Cj^)(CyS+ky)+3 

A = s^+(f2^-3)+(I s+k 
nn ' n y y z 


If only one elastic mode is retained in the model Eq. (4.67) reduces to, 


s^-3fi 

y 

3(1+S7y) 


21 


(n) 


1 


where, A and A, „ 
YY 1,1 


\y 
■(7 s+k ) 

^ y y M-il z 

are defined in Eq. (4.67) 


(CyS+ky)C^^^-2I^^^£/Jy) 

-d+c^ ) (7 s+k^)cf^^+2I,^^^il 
1 y y 2 j" 


1,1 


If the retained mode is an axi-symmetric mode, i.e., 
can be simplified to, ^ 

s®+a,s®+a»s‘*+a„s^+oc, s^+a-s+a, = 0 
1 2 3 4 5 6 

where, a, = (1+c. )c 

■ y 

a„ = (l+c,)k +fi?-3J2 
2 ' 1“^ y 1 y 


= Cy{(l-+c^)(n^-3)-3c^(l+f2^}+4I^^^ 

a, = k {3(l-c )+(l+c,)(n?-3)}+3n^(l-J2 ) - 9 + 

*y *^y xy 


0, Eq* (4.68) 

(4.69) 


41 


( 1 ) 


M, J 
1 y 


C - c 

1 y 
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41 


( 1 )' 


(c,y3) 

a^ - 3(S^J-3)<l-c^ny)ky-9J2y(nJ-3) 

The necessary conditions for stability are given by, 
c. > 0 ; k > ^ 


y (l+Cj^) 
41 


( 1 ) 


nJ-3 > {3c^(14fiy) - } /(1+Cj^) 


y y 


k >{ 9 - 3 fiJ(l-n )- ...... }/{ 3 (x-c,J^^)+(l+c,)(fi?“ 3 )) 

y 1 ' y' 1 y 11 


41 


( 1 ) 


nf-3 >•- 




(Ciky+3) /3(l-CiS2y)Cy 


If, on the otherhand, the retained mode is a mode with one nodal 
for which, l|l) =» 0 , Eq. ( 4 . 68 ) simplifies to, 

s®-l^is®4<!t«s‘*+a_s^+ay s^+a-s+a. ■ 0 
1 2 3 4 5 6 

a./Z . l+.,+C<l-l) 

^ j «L 


a» = k (a. /c )-H 2 ?- 3 n 
2 y' 1 ' y*^ 1 y 

a^/Cy = (S2^-3)(l+c^)-3(l-c^ny)+3C^^’^^ 
» ky(a3/Iy)+3f2j(l-ny)-9 
a^/Cy = 3(J^2-3)(l-flyC^)-9J2yC^^»^^ 


(1-n ) 

N y' 


The necessary conditions for stability are given by, 


c > 0 

y 


k > - (S 2 ,^- 3 f 2 )/(l+c, + C '"’’"0 
y 1 y" 1 


( 4 . 70 ) 
diameter 

( 4 . 71 ) 
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nJ-3 > !3(i-Oj^ny)-3o^^’^^i-ny)Ki+Oj^) 
ky >{9-3S2j<l-t!y)}/(a3/Oy) 
nJ-3 > 3nyC^^’^Va-nyOj) 

ky > 9ny(fJj-3)7(05/Cy) (4.72) 

Wtien more than one elastic mode is retained in the model, a formal expansion 
of the characteristic determinant in Eq. (4.67) is algebraically complicated. 
Hence, we employ a digital algorithm to evaluate the characteristic equation 
when all the system parameters are known. For this purpose Eqs. (4.56), 
(4.57) and (4.58) are re-written in the state variable form as, 

X - Ax (^-73) 

»h«e, . [g Y ___ g, j., ___ g, jT 


0 : I 

h - - 

p ; Q 


0 » null matrix of size (n+2)x(n+2) 

1 ■ identity matrix of size (n+2)x(nH-2) 


y 

-3(l+^2y) 


-{3+k (1+c )} 
k 

y 3 z 


-k . . 

y z 

... 

-(nJ-3)-kyC<^’^^, 


k c 
y 3 z 


-k 

y z 
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• • f 


0 


Q 


. 0 


-d+C3^)Cy (l+c^)CyC^^>-2l(^^ (A/Jy) 

z C -c 

y 3 z y z 


• • • 


» • • 




c C-C 
y 3 z 


(n) 


-c C 
y z 


=3 ■ W 


(n+2)x(n+2) 


One can easily find the characteristic roots of the system of Eqs. (4.73) 
using a digital computer . The result of such an analysis is pre- 
sented in Fig. 4.3. This figure shows the root loci of the lowest frequency 
mode of the system. It can be noted from the figure that the influence 
of the axi-syimetric modes, e.g. (1,0), (2,0) etc., on the other system 
modes is very weak. Since, the coupling between the axi-symmetric modes 
and the rigid body modes is very weak a negligible amount of damping is 
imparted into the axi-syrametric elastic modes. Hence, the characteristic 
roots corresponding to the axi-symmetric nodes lie very close to the 
imaginary axis. We also note from Fig. 4.3 that with an Increase in 
spring stifness, k, the characteristic toots corresponding to the lowest 
frequency modes move toward the imaginary axis. 


B. The case of the nodal diameter along the roll axis (body z axis). 

For this position of the nodal diameter, Pq 0 in Eq. (B.4) in Appendix -B 
Hence , 

C' 


and 


•<“>=. 0 


,(n) 


Vo, P 7* 1 


(l+Cj^p/2 


A study" of Eqs. (4 i 54) -('4. 5 9^) reveals that for the present ^ case*, * the* roil (({)), 
yaw (ijj) , 5 motion and elastic modes with only one nodal diameter form an inde- 
pendent set from that of the pitch (0), y motion and the axi-symmetric elastic 
modes. Hence, we arr.1;v*e at the following two independent characteristic 
equations for the system. 
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Roll-yaw -6 


s "*0. 


(1-S2 )s 


•(1-S2jj) 


s^-4Q_ 


A(l-J2^) 


-(CjjS+k) 


^66 

-(7^s+^)03C<1) 


-(c 3+k )c„C, 


(n) 


z z''^3 y 


Pitch - Y - axi-eyininetric modes 


• * • 


-(l+C2)(rjjS+kjj)C 


1,1 


( 1 ) 


n,l 


• • • 


• • » 


• • • 


(4.75) 


S^-3Qy 

3(l+0y) 


2(I^^Vm^A)s 


•(CyS+ky) 


■2(i['^^Vjy)S 


Ayy 2(I<1U/Jy)s 


s^+nJ-3 


• • • 


. . . 0 


0 9 0 


s^+S7^-3 

n 


(4.76) 


where, ■ s^+(l+C2) (c^s+k^)+4 

A “ (c s+k (mT^n) 

ran z z y 

A = s^+Q^- 3+(7 s+k )C^""^ 
nn n ' z z y 

Ayy = s^+(l+c^)(c^s+ky)+3 

J 

^3 "" M 

n 
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Since a formal expansion of the characteristic determinants in Eqs. (4.75) 
and (4.76) is involved, the numerical algorithm used in the previous cases 
is used to evaluate the characteristic equation and the characteristic roots. 
For this purpose we re-write Eqs. (4. 54) -(4. 59) in the following state variable 
form: 


Xi 

X2 


1,X], 


(4.77) 

(4.78) 


where, " [ 0 Y Eq,! ® 


0,2 


9' Y’ 
o,n ' 


X2 “ [ ^ ®1,2 


e. 

l,n ^ 


r 


0,1 

6' 


o,n- 


1,1 


.e 


l,n 


] 


0, n 

1, n 


non-dimensional modal amplitude of the nth axi-symmetric mode, 
non-dimensional amplitude of modes with one nodal diameter. 


The characteristic roots of the system of Eqs. (4.77) and (4.78) can be 
found by using digital computer algorithm The root loci of the lowest 

frequency mode of the system for the roll-yaw degree of freedom is presented 
in Fig. 4.4. It can be observed that this plot is very similar to the plot 
presented in Fig. 6 of Chapters. *hus, we see that the behavior of a free 
shallow spherical shell is very close to that of a flat circular plate. It 
is also observed that in the pitch degree of freedom, the axi-symmetric modes 
have negligible influent, e on the system characteristic roots correspondi.’g 
to the lowest frequency modes. 
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Fig. 4.1b ; Shallow Spherical Shell with dumbbell 
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Fig, 4.4: Root Locus of Lowest Frequency Mode - Shallow Spherical Shell with Dumbbell 


5A. CONCLUSIONS 


1. Space structures with their maximum moment of inertia axis along 
the local vertical are gravitationally unstable. Such structures 
can be gravitationally stabilized by attaching a light weight 
dumbbell with heavy tip masses. 

2. In the case of the flexible beam, the dumbbell excites only the 
anti-symmetric elastic modes of the beam. 

3. In the case of the flat circular plate, the stabilizing dumbbell 
excites only those elastic modes with one nodal diameter. The 
small concentrated masses added for the roll-yaw stabilization 

of the plate result in only two possible positions for the nodal 
diameters. 

4. In the case of the shallow spherical shell, the axi-syrametric modes 
are coupled to the pitch motion. In addition, the stabilizing dumbbell 
excites those elastic modes with one nodal diameter. The small con- 
centraf“d masses added to stabilize roll-yaw result in only two 
possible positions for the nodal diameter. 

5. The shallow shell undergoes a small static deformation under the 
action of gravity and centrifugal forces. 

6. In all the cases above, as the torsional spring stiffness is increased, 
the characteristic roots of the least damped mode move toward the 
imaginary axis. The stability of the truncated model deteriorates 
with the increase in the number of elastic modes retained in the model. 

1 , To damp the motion of the system in all its modes, especially the 
low frequency modes, the use of active dampers (control systems) 
are needed. Moreover, it is thought that with the use of the passive 
gimballed dumbbell stabilization- device together with active controllers, 
both the peak forces as well as fuel consumption could be significantly 
reduced. Such a study would represent a logical extension of the present 
work. 
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5B. RECOMMENDATIONS 

Possible future contributions to this work could consist of; a stvdy of 
the effects of higher order gravity-gradient potential terms and of solar 
radiation pressure for very large structures* 


In the case of shell structures it was shown that the gravity and orbital 
motion can excite elastic modes of the structure. *Hence, it is suggested 
to investigate the Influence of orbital eccentricity on the elastic 
motion and find possible sources of resonance. 

At the operational altitudes of the future missions involving large 
structures » the principal enviromental disturbance is that due to solar 
radiation. Hence, it is proposed to develop accurate models for solar 
radiation pressure. Thermal gradients are also Induced due to solar 
radiation which can excite some of the lower elastic modes of the struc- 
ture. An investigation of the possible thermally induced structural 
oscillations is recommended. 
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APPENDIX - A 


Kinematic Relations 


Fig (A.l) shows the sequence of Euler angle rotations adopted in going 
from orbit frame, Tq(X 5 ,Yjj,Zq) to the body frame, t^j(X,Y,Z). 



pitch 

roll 

yaw 


Fig. A.l: Euler angle rotations 


The body frame, T|j(X,Y,Z) and the orbit frame, (Xo»Yo»Z^) are related 
to each other by the following transformation relation, 


where , 


b o 



C(j)C0 

s(j)ct|Hn:(})s0st|) 

s(j)st|;-c(j)s0cif^ 

T - 

“S(j)c0 

c(j)cii^-s(i>s9sij^ 

c^stp+stpsQcijj 

X 

S0 

-C0St|J 

cOcip 


T^e,body angular velocity components (£u ,0) ,0) ) and Eular angular rates 
(9,^,i|;) are related as follows, x y z 


to 


X 

to 


y 

to 


z 


0s<|)44*c{J>c0--tOj, (s<j)ct|H"C<|)s9si|^) 

• • 

9c<j)-t|JS({)c9-(Oj, (ctf)cijj-s({)s 9si|j) 
t|^s0+<j)-t(o^c9st|J 
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The transformation between the co-ordinate freme, and the local 
cylindrical co-ordinate frame is as shown below. 



Fig. A.2 




local cylindrical frame 



0 0 
c3 -s3 
s3 c3 


(A. 2) 
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Fig. (A. 2) shows the sequence of Y and 6 rotations adoptijd in arriving 
from body frame, T^j, to the dumbbell frame, T^. 



Fig. A. 3 

s 

The body frame, and the dumbbell frame, T^j, are related to each 
other by the following transformation: 


Tb - ^2 td 



where. 





cYc6 

-cys6 

sy 


s6 

c6 

0 

2 

-syc6 

sys6 

cy 


(A. 3) 
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APPENDIX B 


Frequencies and Mode Shapes of a Shallow Spherical Shell*; 

The results presented here are based on the following assumptions. 

• the shell has constant thickness 

• the vibration of the shell is primarily in the transverse direction 
i.e. parallel to shell's symmetry axis. The effect of longitudinal 
inertia is negligible in comparison with transverse inertia. 


The natural frequencies of a thin, shallow spherical shell with completely 
free edge can be obtained from the roots (Xi ,X2 n****) of tran~ 
scendental equations, 


0) >w 
n <» 


W <0) 
n « 

where, 


x‘ 

1?^ Rp(A) 


+ 


ipW 1 


1-V 


(p«0,l) 


- 1 (p-2,3,.. . .) 




Up(n) 

Vp(n) 

I 

n 00 ^ 


(p“2 , 3 p . . . ) 


(j“ 1,2,...) 


(B.l) 


(B.2) 


(B.3) 


* natural frequency 
I » base radius of the shell 
h » wall thickness of the shell 

p = mass per unit volume of the material of the shell 
D = flexural rigidity, EhVl2(l-V^) 

V ^ Poisson's ratio 
(D„ » (C/hpR^)^ 


C = longitudinal stiffness factor. Eh 

R =» radius of curvature of the middle surface of the shell. 


J (X),I (X) = Bessel function and modified Bessel function of the 

^ ^ first kind and order, p. Primes denote the derivative 

with respect to the argument. 

SpW = 4p2(p^-XHl-v){X[Jp(A)r(X)-J^(X)lp(X)j 

+ (p+l)(l-v)[I^(X)- t i^(j)][j'(X)- £ Jp(X)]} 

* ^ 1 

^Extracted from Johnson, M.W. and Reissner, E. , "On Transevarse VibrUw 
of Shallow Spherical Shells," Quart. Appl. Math., Vol. 15, No. 4, Jan 1958, 
pp. 367-384. 



R (X) - {a-v)[xj' (X)-p=*j (X)]-fX'‘j„(X)H(i-v)?Mxi’(x) 

p r P r r 

-Ip (X) 1 - X ' I ’ (X)} -{ (1-v )p 2 [X J V(X) "Jp (X) ]+X ® J ’ CX) } 

x{(l-v)[Xl’(X)-p=^Ip(X)]- X^ipCX)} 

n - (i -/ -X , n >o) 

U (n) •• 2(l-v)p’(p^-I){2^/^pn(ber’nbei n-ber nbei'n) 

P , 2 P P p P 

+ 2^(1-V)p (p+1) [ ^ (ber^n+bei^n) - ^ (ber^n 

+ bei^n) ’+(ber'ri)^+(bei’ri)^] ) 

P P P 

Vp(n) * [(l-V^)p^(p^-X)-n'*]2%(ber^nbeipTi~berpnbei^n) 

+ 2^(l-V)nHp^[ “ (ber^n+bei^n)]’-<ber'n)^ 

n p p p 

- (bei’n)^} 

3/2 

ber T) “ real part of J_(i ' n) 

^ 3/2 

beipfi “ imaginary part of Jp(i n) 

k" » a'*/DR^ 

The elastic mode shapes of the shallow spherical shell are given by the 
expression, 




(n) 


»r 


j.p 




(B.4) 


where , 


P 

j 



3 

^0 


= Number of nodal diameters in the n*^^ mode 
= Number of nodal circles in the n*^^ mode 
» Roots obtained form eqns. (E.l), (B.2) and 
“ Non-dimensionalized radial distance, 0£C£l 
=* angle defined in fig. B.l 
= arbitrary phase angle 


(B.3) 
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C. ,D. *• constants determined from Eqns. (6,l)-(6.4) of 

j>P J»P Johnson and Reissner (1958) 


Aj p " arbitrary constant 


Thus, from Eq. (G.13) we can observe that the mode shapes of shallow 
spherical shells are characterized by a nodal pattern consisting of 
a set of concentric nodal circles (j) centered on symmetry axis and 
a set of nodal diameters (0) , The modes with mode shapes having no 
nodal diameter (0"O) are called ’'axi-symmetric" modes. Thus, axi- 
symmetric mode shapes are dependent only on 5 , 


The expressions for the constants are given by. 


.(n). 

'y 3? 


Similarly, 
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,(n). 
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■ 0 (p)*l) 
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(p-i) 


- 0 (p^i) 
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APPENDIX C 


Evaluation of Integrals I^^^ 1^^^ and I^^^; 

x_ Xn r 


Appendix B we have, 


FromEq. (B.13) of 


)P+4 


" ^.P °j,p ^ 

Also, from the equation of a shallow spherical shell 


X (?) - 

c 

By definition, 

(j) 


2R 


d-c") 


V ^ /^x^(?)Wj^o(?)?dC 


(c.l) 


(C.2) 


i.e. 

^Xq^ "‘IF’ 

Using the values of integrals given in Table “Cl and noting that 
Dj^Q ■ ^q)/I^(Xj ^q) we obtain, 


(0.5) 


Now consider. 


(j) 


X.. 


/' X (?)W, o(?)?d? 
0 ^ 


i.e. 


(j) 


X. 




Hence using the values of integrals listed in Table Cl we get, 

i6 D 
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X, 


J2,2 

2r" ^,2 


^ URD^** 2 ^J.2 ^~x\ 2 

J » ^ j 


4 . 1 
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Similarly, consider 

I«). j-1 

r o 3 , j 




J »>^ 


Using the values of the Integrals given in Tabled, we arrive at, 
J » 3 


J » 


(C.5) 
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Table Cl: Some Useful definite integrals. 


0 

0 ^ 

/^cj„ac)dc 
0 ^ 

/^c^j.aodc 

0 

/'c2j„0\C)d? 
0 1 

/^ClQ(AC)di; 

(X?)d? 

0 ^ 

;^?i (A?)d^ 

0 ^ 

/^C^i.C^OdC 

0 

0 

(XOdC 

0 


Jl(X)A 

[X^ (X)+2X Jq (X) -4 (X) ] /X ® 

"[XJ3^(X)+2Jq(X)-2]/X'^ 

[8 Jj^(X) -4X Jq (X) -X^ (X) ] /X^ 

[8-8Jj3 (X) -4X (X) -X"^ (X) ] /X 

Ij^(X)/X 

[X=I^(X)-2XIq(X)+4Ij^(X)]/X^ 

[Xl^(X)-2Io(X)+2]/X" 

[X2Ij^(X)-4XXq(X)+8Ij(X)1/X' 

1 

n+1 

[\\ (X) ~4Xlj^ (X)+8 Iq (X) -8] /X ® 
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